Magnetic Coupling of a Rotating Black Hole 
with Advection-Dominated Accretion Flows 



Yong-Chun Ye, Ding-Xiong Wang 1 

Department of Physics, Huazhong University of Science and Technology, Wuhan, 
\ 430074, People's Republic of China 

^ ; and 
^ : Ren-Yi Ma 



Shanghai Astronomical Observatory, Chinese Academy of Sciences, Shanghai, 
200030, People 's Republic of China 



> 

qq . Abstract 

A model of magnetic coupling (MC) of a rotating black hole (BH) with advection- 
dominated accretion flow (MCADAF) is proposed. It turns out that MCADAF 
providers a natural explanation for the transition radius between ADAF and SSD, 
and could be used to interpret the highest luminosity of GX 339-4 in hard-state. A 
very steep emissivity index can be produced in the innermost part of the MCADAF, 
which is consistent with the recent XMM-Newton observations of the nearby bright 
Seyfert 1 galaxy MCG-6-30-15 and with two X-ray binaries (XRBs): XTE J1655-500 
and GX 339-4. In addition, we estimate the BH spins in Seyfert 1 galaxy MCG-6- 
30-15 and in the two XRBs based on this model. 

Key words: accretion, accretion disk, black hole physics, magnetic field, 
hydrodynamics — 97.60.Lf, 04.70.-s, 95.85.Sz, 95.30.Lz 



1 INTRODUCTION 



More than three decades ago, the standard thin disk was proposed by Shakura 
& Sunyayev (1973, hereafter SSD), which is regarded as the milestone of black 
hole (BH) accretion theory. However, the effective radiation temperature is 
very low in SSD, being obviously not consistent with the observations of the 
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hard X-ray and gamma-ray radiation in X-ray binaries (XRBs). In order to 
overcome this shortcoming, some authors proposed the cooling-dominated hot 
disk (Shapiro, Lightman & Eardley 1976, hereafter SLE). A vital disadvantage 
with SLE lies in thermal instability. Not long ago, advection-dominated ac- 
cretion flow (hereafter ADAF) was suggested by some authors (Narayan & Yi 
1994, hereafter NY94; 1995a; 1995b; Abramowicz 1995; Esin 1997; Manmoto 
1997; Gammie and Popham 1998; Popham and Gammie 1998; Yuan 1999). 
The ADAF model succeeds in explaining the observed spectra of some sources, 
such as Sgr A* (Narayan 1995c), NGC4258 (Lasota 1996) and GRO J1655-40 
(Lasota 1998). Recently, Yuan (2001 hereafter Y01; 2003) extended the ADAF 
model and put forward a new two-temperature hot branch of equilibrium so- 
lutions, i.e., the luminous hot accretion flow (hereafter LHAF). It is argued 
in Y01 that LHAF is much more luminous than ADAF as the entropy of the 
accreting mater supply the radiation as well as the viscous dissipation. 

Yuan (2006, hereafter Y06) pointed out that the hard states at their high- 
est luminosities of some XRBs could not be explained by ADAF, while LHAF 
can produce the highest luminosity in the hard state of XTE J1550-564. How- 
ever, as argued in Y06 that neither ADAF nor LHAF could be used to explain 
the highest luminosity of GX 339-4 (Done & Gierlinski 2003; Zdziarski et al. 
2004), which is a great challenge to LHAF. Probably, an efficient way of in- 
creasing the luminosity in hard-state is the magnetic coupling (MC) process, 
which is an important mechanism for extracting energy from a rotating BH 
(Blandford 1999; Li & Paczynski 2000; Li 2002a, hereafter L02; Wang et al. 
2002, 2003, hereafter W02, W03, respectively). However, the MC has been for- 
mulated in the context of a relativistic thin disk (henceforth MCTHIN, Page 
& Thorne 1974; L02; W03), in which both the gravitational energy of the ac- 
creting matter and the rotational energy transferred into the disk are assumed 
to be radiated away in black-body spectrum. Thus MCTHIN is obviously not 
consistent with observed hard power-law X-ray spectra. 

To solve the above puzzle in LHAF we intend to incorporate the MC with 
ADAF (henceforth MCADAF). It turns out that the steep emissivity index in 
Seyfert 1 galaxy MCG-6-30-15 and the two XRBs: XTE J1655-500 and GX 
339-4 could be interpreted based on MCADAF, and our fittings are consistent 
with the low BH spins in these sources, and that the observation of the highest 
hard-state luminosity of GX 339-4 could be explained. Furthermore, a tran- 
sition radius between ADAF and SSD could be given naturally in MCADAF 
with the hard X-rays spectra produced in the MC region. 

This paper is organized as follows. In § 2 MCADAF is outlined, and the 
expressions of the MC power and torque in MCADAF are derived. In § 3, 
the MC effects on disk radiation and the emissivity index are discussed. It is 
found that the emissivity index is consistent with the recent XMM-Newton 
observations of the nearby bright Seyfert 1 galaxy MCG-6-30-15 and the two 
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Fig. 1. The poloidal magnetic field for MCADAF, in which ADAF and SSD are 
located inside and outside of r out , respectively. 

XRBs: XTE J1655-500 and GX 339-4 for the low BH spin. Finally, in § 4, 
we summarize our main results. Throughout this paper the geometric units 
G — c — 1 are used. 



2 DESCRIPTION OF MCADAF 



In the previous MC model the disk is regarded as a thin disk, in which 
the rotational energy transferred from the BH to the inner disk is assumed 
to radiate away as the black body spectrum. In MCADAF, however, this as- 
sumption is changed to overcome the shortcomings of MCTHIN. It is assumed 
that only a fraction of energy transferred in the inner disk is radiated, and the 
rest energy deposited in the accreting matter as entropy is delivered to the 
BH as ADAF. Thus we have the magnetic field configuration of MCADAF as 
shown in Figure 1, in which ADAF is located in the MC region with closed 
field lines penetrating the accretion flow, and SSD is located outside the MC 
region. This configuration is significantly different from that given in W02 for 
MCTHIN, and a transition radius between ADAF and SSD is given naturally 
in MCADAF. 
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In order to facilitate the discussion on MCADAF in an analytic way, we 
make the following assumptions. 



(i) ADAF is both stable and perfectly conducting, and the closed magnetic 
field lines are frozen in the accretion flow. 

(ii) The magnetic field is assumed to be weak, and the instabilities of the 
magnetic field are ignored. The magnetosphere is stationary, axisymmetric 
and force- free outside the BH and ADAF. 

Following Armitage & Natarajan (1999), the total pressure P in the inner 
part of MCADAF is given by 

P=^(5 + 2e') 1/2 (M/r 5 )^, (1) 



where the magnetic pressure P mag is related to the total pressure P by 

Pmag = ^/StT = (1 " /3 S )P. (2) 

In equations (1) and (2) a s is the viscosity parameter, (3$ is the ratio of gas 
pressure to total pressure, and e' is defined as e' — e/f — jj^h (NY94). The 
parameter / measures the degree to which the flow is advection-dominated, 
and 7 = fzfff is ^ ne ratio of specific heats (Esin 1997). Thus we have 

,<=^. (3) 



Combining equations (1) and (2), we have the expression for the magnetic 
field in the inner part of the MCADAF as follows, 



B D = v / (2v / 2/3)(5 + 2e') 1 / 2 a~ 1 {l - [3 s )M l / 2 Mr- 5 / 2 . (4) 

Narayan et al. (1998) pointed out that the inner edge radius of an ADAF 
is located between r mb and r ms , i.e., 

^'mb — Tin — ^'ms- (5) 

where r mb and r ms are the radii of the innermost bound circular orbit and the 
last stable circular orbit, respectively. The radius r in can be expressed by 

r in = MXi n , 

(6) 

"X.in Xmb ^(Xms Xmb): 
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where Xmb = yJr mb /M and Xms = \jr ms jM. The parameter A is used to 
adjust the position of the inner edge of the ADAF, and A = 0.5 is adopted in 
calculations. Thus equation (4) can be expressed by 

B D » B CC 5/ \ (7) 
B = M- x M xj2 = 1.38 x l^m-^m 1 12 Gauss, (8) 

C = (2v / 2/3) 1 / 2 (5 + 2e') 1/4 a; 1/2 (l - Ps) 1 ^ 2 , (9) 



where £ = r/r in is a dimensionless radial parameter, is the BH mass 

in the units of one solar mass, and rh = M j {L E j (?) is the accretion rate in 
terms of the Eddington luminosity. 

Incorporating equations (3) — (9), we can calculate the magnetic field B D 
in MCADAF versus £ with the given a* and / as shown in Figure 2. In the 
following calculations, rh = 0.01, a s = 0.3, @s — 0.5 are adopted as given by 
Narayan et al. (1998). 

The magnetic field in MCTHIN varying with the radial parameter £ is 
given in W03 as follows, 

B p D = B p H [rJzu D (r ms )]C n - (10) 



Incorporating equation (10) with some relations given in W03, we have 
the curves of the magnetic field B P D versus £ with given a* and n in MCTHIN 
as shown in Figure 3. 

Inspecting Figures 2 and 3, we find that the magnetic field in MCADAF 
is weaker than that in MCTHIN. 

The mapping relation between the angular coordinate on the horizon and 
the radial coordinate in MCADAF can be derived based on the conservation 
of magnetic flux in an analogous way given in W03, and we have 

sin0d0 = -R(a„0df, (H) 



where 



p / n r 1/4 xfn V / i + ^ 4 e- 2 + 2a^~ 

R(a*,4) = — , , (12) 
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Fig. 2. The curves of lg [Bp (Gauss)] in MCADAF versus £ with m BH = 10 for (a) 
/=0.5 and a*= 0.3, 0.5 and 0.7 in solid, dotted and dashed lines, respectively, and 
for (b) a*=0.5 and /= 0.3, 0.5 and 0.7 in solid, dotted and dashed lines, respectively. 

Integrating equation (11) and setting £ = £ in — 1 at 9l — 0.457T, we have 

cos9m — cos9l= J R(a*,£)d£. (13) 



The parameter £ out = r out /r in in equation (13) represents the outer radius 
of the MC region, and we have the curves of £, out versus a* with different values 
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Fig. 3. The curves of magnetic field lg[Bn {Gauss)] in MCTHIN versus £ with 
n = 5/4 and a*= 0.3, 0.5 and 0.7 in solid, dotted and dashed lines, respectively. 

of 9m as shown in Figure 4. 

From Figure 4 we find that the variation of ^ out is insensitive to the value 
of 9m i and the MC region in MCADAF is restricted to a very small range, 
1 < £ < (£out)max = 1-595. This result implies that the transition radius in 
MCADAF is very small. 

Also analogous to W02 we can derive the expression for the MC power and 
torque in MCADAF by using an equivalent circuit in BH magnetosphere (Mac- 
donald & Thorne 1982), and the torque and power exerted on the MCADAF 
between two adjacent magnetic surfaces are given by 



AT MC = zuB H IAl = (A^/2tt) 2 (Q h - Q D )/AZ H , 

ap mc = n F AT = (a^/2tt) 2 n F (n H - n D )/Az H , 



(14) 
(15) 



where 

AZ H = R H Al/(2nw) = 1pA9jw. 



(16) 



Incorporating the expression for angular velocity given in NY94, we ex- 
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Fig. 4. The outer radius of the MC region, £ out , versus a* with 6m= 0, 0.057T and 
0.1 7T in solid, dotted and dashed lines respectively. 

press the angular velocity VL D in MCADAF as follows, 

Q D » [2e'/(5 + 2e')] 1/2 fi*r, V K = * /2 3 r - (17) 



Integrating equations (14) and (15) over the angular coordinate from 6m 
to 6*l = 0.457T, we obtain the total MC power and torque as follows. 



T /T 4a*C 2 xl n (l + xIr?< + 2xX) f (1 -P) sin* OdO 

±MC J-0 = 7—. / ~ T\ \ ■ 2/) ' l 18 J 

1 + q J 2 — (1 — q) sin 9 

9m 



PmcI Pq 



2alC 2 X j n (l + X " 4 a 2 + 2 X ~Sl) P) sin 3 Odd 



/■ 



g) sin 2 ' 



(19) 



where 



^W ^'^'^.y , (20) 
a, (a, + C 3 / 2 X 3 „) 
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Fig. 5. The curves of ^ out (solid line) and £ c versus a* with /=0.3, 0.5, 0.7 in dashed, 
dotted and dot-dashed lines, respectively. 

P = (B 2 ) M 2 ~ m BH m x 1.25 x 10 36 erg ■ s' 1 
T = {Bl) M 3 « m| H m x 5.52 x 10 29 g • cm 2 ■ s~ 2 

The co-rotation radius r c can be determined from the equation (20) by 
setting (3 = 1, which indicates the place at ADAF with = Qjj. Incor- 
porating equations (12), (13) and (20), we have the curves of dimensionless 
co-rotation radius £ c = r c /r in and £ ou t versus a* for the different values of / 
with 9m — as shown in Figure 5. 

Inspecting Figure 5, we find that each curve of £ c decreases monotonically 
with a*, and it intersects with the curves of ^ out and £j n = 1, respectively. 
These intersections imply that the co-rotation radius r c is located at the outer 
and inner boundary of the MC region for the BH spin equal to a° ut and ct* n , 
respectively. 

The MC region is divided into two parts by £ c , i.e., the inner MC region 
for 1 < £ < £ c and the outer MC region for £ c < £ < £ oui with a° ut < a* < a™. 
Thus the energy and angular momentum are transferred by the MC from the 
BH into the outer MC region with < Qjj, while the transfer direction 
reverses for the inner MC region with > Qjj. The correlation of the BH 
spin with the transfer direction is given as shown in Table 1. 
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Table 1 



The correlation of the BH spin with the transfer direction of energy and angular 
momentum 



Parameters 


Transfer direction of energy and angular momentum 


/ 




out 


a% ut < a, < af 


af < a* < 0.998 


K £ < & 


ic < i < (.out 


£ c does not exist 


0.3 


0.279 


0.175 


from disk to BH 


from BH to disk 


from BH to disk 


0.5 


0.228 


0.146 


0.7 


0.196 


0.127 



By using equations (11) — (21), we have the curves of the MC power and 
torque versus a* for different values of / as shown in Figure 6. It is found that 
the MC power increases monotonically with the increasing BH spin, while the 
MC torque increases generally with the BH spin, except the spin approaching 
unity. It is noted that both the MC power and torque are insensitive to the 
values of the parameter /. 



3 STEEP EMISSIVITY INDEX PRODUCED IN MCADAF 



Based on the MC effects on radiation flux we can calculate the emissivity 
index in MCADAF. The energy conservation equation given in NY94 is 

+ ^_ _ 2a s pc 2 s r 2 H dtl 2 ^_ _ f 2a s pc 2 r 2 H f dQ^ 



where Q + is the energy input per unit area due to viscous dissipation, and Q~ 
is the energy loss due to radiative cooling. Considering the electromagnetic 
flux of energy transferred by the MC process, we modify equation (22) as 
follows, 

2apc 2 s r 2 H dQ 2 f 2a s pc 2 s r 2 H dQ 2 

— + HM ° n - Q =f n K ( ^ } ' (23) 



where H MC is the angular momentum flux transferred in MCADAF, which is 
related to T MC by (Li 2002b, W03) 



dT MC /dr 



AtttHmc 

13) sin 3 I 



_ / 4T q»(l 

^ 2-(l-g)sin 2 



c 2 xL(i+x. 



" 4 „2 



(80) 
\drj ' 



(24) 
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(b) 

Fig. 6. The curves of Pmc/Pq (upper panel) and Tmc/Tq (lower panel) versus a* 
with /=0.3, 0.5, 0.7 in solid, dotted, dashed lines, respectively. 



and 86/ dr can be calculated by using equation (12), 



de/dr = {86/ di) (d£/dr) = ~ ^r^j - (25) 
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Incorporating equations (24) and (25), we have 

A(a*,£) , 1 < f < £ out , 



H(a*,$/H ={ 



(26) 



where 

# = (^)M = 77i x 3.15 x 10 20 g- s~ 2 , (27) 
A ( °*' " lTg TTpcsc^-Cl-g)] 72 K • (28) 



Then we obtain the expression for radiation flux in MCADAF as follows, 
Ftotai = Q~ = Wl + (1 - fpW&^Lf. (29) 



Combining equations (24) — (29), we have 



In equation (30) we have i<o = (Bq) c, and the first and second terms at 
RHS represent the MC effect on radiation flux and the contribution due to 
ADAF, respectively. 

Following L02 and W03, the emissivity index is defined as 

a = -d\nF/d\nr. (31) 

By using equations (24) — (31), we have the curves of the emissivity index 
a versus lg (r / M) for the different values of a* and / as shown in Figure 7. As 
shown in the shaded region of Figure 7, we find that the recent XMM-Newton 
observation of the nearby bright Seyfert 1 galaxy MCG-6-30-15 (Wilms et al. 
2001) can be simulated in MCADAF, while the index produced by the ADAF 
is far below the shaded region as shown by the dot-dashed lines in Figure 7. 

Comparing the results in Figure 7 with those obtained in W03, we find 
that the emissivity index in MCADAF corresponds to the spins less than 0.4, 
which are much lower than those obtained in MCTHIN with a* > 0.99. By 
analyzing the X-ray spectrum of the Seyfert 1 galaxy MCG-6-30-15, Dovciak 
et al. (2004) discussed a low-spin (a* = 0.25) best-fitting model for the iron 
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Fig. 7. The emissivity index of the Seyfert 1 galaxy MCG-6-30-15 inferred from the 
observation of XMM-Newton is shown in the shaded region. The emissivity index 
in MCADAF is plotted versus lg(r/M) for a*=0.30 and 0.35 in panels (a) and (b), 
respectively. The curves are plotted for /=0.3, 0.325 and 0.35 in solid, dotted and 
dashed lines, respectively. The emissivity index in ADAF without MC is plotted in 
dot-dashed line. 

line emission in MCG-6-30-15. This result implied that the central BH spin 
may be low. Brennemanl and Reynolds (2006) pointed out that an inner 
emissivity profile requires an inner radius of r inner = 3.2r g (r g = M) to fit the 
iron line emission. These results are consistent with MCADAF but not with 
MCTHIN, since the inner edge of MCADAF is located at r in < 3.2r g for a 
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low BH spin, and the inner emissivity profile can be produced by virtue of the 
energy transferred in the MC process. 

It is shown by the recent observations that the very steep emissivities are 
also found in XTE J1650-500 (Miller et al. 2002; Miniutti et al. 2004) and GX 
339-4 (Miller et al. 2004). The simulations for the steep emissivity in these 
XRBs can be worked out in the same way as MCG-6-30-15, and the curves of 
the emissivity index a versus lg (r/M) for the different values of a* and / are 
shown in Figure 8, from which we find that the BH spins of the two XRBs are 
also estimated as the low values in fitting the observed very steep emissivity 
indexes. 



4 DISCUSSION 



In this paper, we propose a toy model for advection-dominated accretion 
flows around a rotating BH with the global magnetic field. And the expressions 
for the extracting power and torque in the MC process and the MC effects 
on radiation flux are derived. By simulating the emissivity index in Seyfert 1 
galaxy MCG-6-30-15 and the two XRBs: XTE J1655-500 and GX 339-4, we 
find that the steep emissivity index could be interpreted based on MCADAF, 
and our fittings are consistent with the low BH spins in these sources. 

Comparing with MCTHIN and ADAF (or LHAF), the advantage of MCADAF 
lies in the following aspects. 

(1) The electrons of high-temperature in the optically thin plasma are 
contained in MCADAF, and they are cooled via synchrotron, bremsstrahlung 
and inverse Compton processes, being responsible for producing hard X-rays 
spectra in a natural way. 

(2) MCADAF providers a natural explanation for transition radius be- 
tween ADAF and SSD as shown in Figure 1. 

(3) MCADAF could be used to interpret the highest luminosity of GX 
339-4 in hard-state. 

Following W02, we have the energy extracting efficiency in MCADAF as 
follows: 

Vmcadaf Vadaf Vmc 

(32) 



14 



1 



0.63 0.66 0.69 0.72 

lg(r/M) 

(a) 



a*=0 . 35 




/ 


GX 339-4 


// 




// 




// 




// 




*/ 


















1 1 





0.62 0.65 0.68 0.71 

lg(r/M) 

(b) 

Fig. 8. The emissivity indexes of XTE J1650-500 (upper panel) and GX 339-4 (lower 
panel) inferred from the observation of XMM-Newton are shown in the shaded 
region. The emissivity indexes produced in MCADAF versus lg(r/M) are plotted 
in upper panel in solid, dotted and dashed lines for a*=0.30 with /=0.25, 0.275 and 
0.3, respectively. Those are plotted in lower panel in solid, dotted and dashed lines 
for a*=0.35 with /=0.18, 0.2 and 0.22, respectively. The emissivity index produced 
in ADAF without MC is plotted in dot-dashed line. 

where 



Vadaf ~ ^ini 

Vmc = Pmc/m, 

E in = (l- 2 X ~n + a^xl. 



(33) 



1 - 3x" 2 + 2a,xr„ 3 ) 
15 
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Fig. 9. The curves of rj ADAF (solid line) and rj MCADAF (dashed line) versus a* with 
/=0.5; (b) the ratio Vmcadaf/Vadaf versus a* with /=0.5. 

By using equations (32) and (33), we have the curves of r] MCADAF and 
Vadaf versus a *j an d the ratio of the former to the latter versus a* as shown 
in Figure 9. 

From Figure 9, we find that the energy extracting efficiency in MCADAF 
is greater than that in ADAF for a* > a* n , and the former could be 1.88 
times at most than the latter for a* = 0.9. This result implies that the lu- 
minosity could be augmented significantly in MCADAF. As argued in Y06 
the maximum luminosity in LHAF is L^haf ~ 16%Lg, then the maximum 
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luminosity in MCADAF could be Lmcadaf ~ 30%Le, which could be used 
to explain the observation of the highest luminosity of GX 339-4 in hard-state 
with 25 - 30%L E (Done & Gierlinski 2003; Zdziarski et al. 2004). 

In this paper MCADAF is worked out by combining the MC process 
with the self-similar solutions given in NY94. However, there exists an obvious 
inconsistency in this model. The self-similar solution is given based on pesudo- 
Newtonian potential in NY94, while the MC is formulated in the context of 
general relativity. Not long ago, some authors (Gammie and Popham 1998; 
Popham and Gammie 1998) obtained the precise solutions in ADAF in the 
context of general relativity, and these works provider a possibility to improve 
MCADAF in our future work. 
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